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H ■ Abstract 

The distribution of the return intervals r between price volatilities above a threshold height q 



for financial records has been approximated by a scaling behavior. To explore how accurate is 
the scaling and therefore understand the underlined non-linear mechanism, we investigate intraday 
datasets of 500 stocks which consist of the Standard & Poor's 500 index. We show that the 



> 

C^) . cumulative distribution of return intervals has systematic deviations from scaling. We support 

this finding by studying the m-th moment /i m = ((r/ ' (r))' m ) l ' m , which show a certain trend with 

t^! ■ 

the mean interval (r). We generate surrogate records using the Schreiber method, and find that 
their cumulative distributions almost collapse to a single curve and moments are almost constant 
for most range of (r). Those substantial differences suggest that non- linear correlations in the 
' original volatility sequence account for the deviations from a single scaling law. We also find that 

the original and surrogate records exhibit slight tendencies for short and long (r), due to the 
discreteness and finite size effects of the records respectively. To avoid as possible those effects for 
testing the multiscaling behavior, we investigate the moments in the range 10 < (r) < 100, and 
find the exponent a from the power law fitting \x m ~ (r) a has a narrow distribution around 
which depend on m for the 500 stocks. The distribution of a for the surrogate records are very 
narrow and centered around a = 0. This suggests that the return interval distribution exhibit 
multiscaling behavior due to the non-linear correlations in the original volatility. 
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I. INTRODUCTION 



The price dynamics of financial markets has long been a focus of economics and econo- 



physics research 

n a a a q, b a a. 

Studying the volatility time series is not only 
crucial for revealing the underlined mechanism of financial markets dynamics, but also 
useful for traders. For example, it helps traders to estimate the risk and optimize the 



portfolio 8 

baa 




The volatility series is known to be long-term power-law correlated 



17, 



18 



19 



2(J, |2l|. To better understand these correlations 



and characterize temporal scaling features in volatilities, recently Yamasaki et al. 
Wang et al. [23 



22j and 



24[ studied the statistics of return intervals r between volatilities that are 



above a given threshold q, which is an alternative way to analyze long-term correlated time 
series (see Ref 25j and references therein) . They find that scaling and memory in the return 



intervals of daily and intraday financial records are similar to that found in the climate and 



earthquake data 



26 



28]. 



Studies of financial records show that the scaling in the return intervals distribution can 



be well approximated by a scaling function 



22 



23, 



24|. However, financial time series 



are known to show complex behavior and are not of uniscaling nature [29| and non-linear 



features 1 301 . Recent studies 



31 



32 



33| of stock markets show that the distribution of 



activity measure such as the intertrade time has multiscaling behavior. Thus, a detailed 
analysis of the scaling properties of the volatility return intervals is of interest. It might 
improve our understanding of the return intervals statistics and shed light on the underlined 
complex mechanism of the volatility. Our analysis suggests that for all Standard & Poor's 
index constituents, the cumulative distributions of the return intervals depart slightly but 
systematically from a single scaling law. We also find that the moments fi m = ((r/ (r))' m ) 1 / m 
are consistent with the deviations from scaling. However, using the corresponding surrogate 



records 



34 



35 



361 ] which remove the non-linearities, \i m almost does not depend on (r) and 



no deviation from scaling occur. Therefore, our results suggest that non-linear correlations 
in the volatility account for the deviations from a scaling law. 

The paper is organized as follows: In section II we introduce the database and define 
the volatility. In section III we discuss the scaling and investigate the deviations from 
scaling in the cumulative distributions of the return intervals. We also describe the stretched 
exponential form suggested for the distribution and the generation of the surrogate records. 
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Section IV deals with the moments of the return intervals. We quantify the deviation 
from the scaling that exhibits multiscaling behavior. We simulate the return intervals with 
different sizes and show the finite size effect for long (r). We also study the discreteness 
effect for short (r) and explore the relation between the moment and its order. In Section 
V we present a discussion. 



II. DATABASE 

We analyze the Trades And Quotes (TAQ) database from New York Stock Exchange 
(NYSE), which records every trade for all the securities in United States stock markets. 
The duration is from Jan 1, 2001 to Dec 31, 2002, which has a total of 500 trading days. We 



study all 500 companies which consist of the Standard & Poor's 500 index (S&P 500) 37], 
the benchmark for American stock markets. The volatility is defined the same as in 23]. 
First we take the absolute value for the logarithmic price change, then remove the intraday 
U-shape pattern, and finally normalize it with the standard deviation. Here the price is the 
closest tick to a minute mark. Thus the sampling time is 1 minute and a trading day usually 
has 391 points after removing the market closing hours. For each stock, the size of dataset 
is about 200,000 records. 



III. SCALING IN RETURN INTERVALS 



The probability density function (PDF) for the return intervals r of the financial volatil 
ities is well-approximated by the following form, 

1 



. 3 



(r) 



/(r/<r», 



(1) 



as analyzed by Yamasaki et al. [22j and Wang et al. [23|, |24|. Here (•) stands for the average 
over the dataset and t/(t) depends on the threshold q. It was s ugg ested that the scaling 



function can be approximated by a stretched exponential [2^ 



23 



24) 



/(*) 



ce 



-(ax)~i 



(2) 



for financial records, which is consistent with other long-term correlated records 



25 



26 



27j, 



281 ]. Here a and c are fitting parameters and 7 is the exponent characterizing the long-term 
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correlation 23, |26|, |27|, [28[. From the normalization of PDF follows 38], 

1= / P q {r)dT= / -i-ce-^/^dr. (3) 
Jo Jo {V 

From the definition of (r) follows, 



(r) 



f*00 /*OC 1 

/ r-P q (r)dr= / r • — ce- (aT/<r))7 c/r. (4) 
Jo Jo {V 

Thus, using Eqs. ([3]) and (jlj), the parameters a and c can be expressed by 7, 

a = r(2/ 7 )/r(i/ 7 ), 

c = 7 a/r(i/ T ) = 7r(2/ 7 )/r(i/ 7 ) 2 . (5) 

Here T(a) = J °° t a ~ 1 e~ t dt is the Gamma function. Thus, if the stretched exponential dis- 
tribution is valid for the scaled interval t/(t), it is completely determined by 7. For 7 = 1, 
the record has no long-term correlations and the return interval distribution indeed follows 
an exponential distribution, represented by a Poissonian statistics, as expected. 

Though the scaling in the return intervals distribution is a good approximation, we find 
slight deviations that as shown below are attributed to non-linear features. To explicitly 
explore the quality of the scaling in return interval distributions, we study all S&P 500 con- 
stituents and show the results of four representative stocks, Citigroup (C), General Electric 
(GE), Coca Cola (KO) and Exxon Mobil (XOM). All other stocks studied here usually show 
similar features. First, we examine the cumulative distribution of the scaled intervals. 

POO POO 

D(t/(t)) = / P q (r)dr = / f(x)dx. (6) 

Jt Jt/{t) 

If the scaling function f(r/ (r)) is valid, the cumulative distributions should also collapse 
to a single curve. Otherwise, the cumulative distributions, which integrate deviations, may 
show clearer deviations from scaling. Indeed, in Fig. [1] we show cumulative distributions 
for three thresholds q = 2, 4 and 6. Note that the volatility is normalized by its standard 
deviation, the threshold q is in units of standard deviations and therefore q = 6 is a quite 
large volatility. It is clearly seen that those distributions are close to each other but do 
not collapse to a single curve. More important, they show apparent deviations from the 
scaling, which are systematic with the threshold. For small scaled intervals (t/(t) < 1), 
the cumulative distribution decreases with q, while for large scaled intervals (r/(r) > 1), 



it increases with q 



In other words, the scaled interval prefers to be larger for higher 



threshold. This systematic trend suggests multiscaling in the return intervals, which might 
be related to the non-linear correlations in the volatility. 

To better understand the systematic trends and test if it is not due to finite size effect or 
discreteness of minutes, we also measure the cumulative distribution of return intervals for 
surrogate records of volatilities using the Schreiber method 34], |35|, |36( where non-linearities 
are removed. For a given time series, we store the power spectrum and randomly shuffle the 
sequence, then we apply the following iterations. Each iteration consists of two consecutive 
steps: 

(i) We perform the Fourier transform of the shuffled series, replace its power spectrum 
with the original one, then take the inverse Fourier transform to achieve a series. This 
step enforces the desired power spectrum to the series, while the distribution of volatilities 
usually is modified. 

(ii) By ranking, we exchange the values of the resulting series from step (i) with that of 
the original record. The largest value in the resulting series is replaced by the largest one 
in the original series, the second largest value is replaced by the second largest one, and so 
on. This step restores the original distribution but now the power spectrum is changed. 

To achieve the convergence to the desired power spectrum and distribution, we repeat 
these two steps 30 times. By this way, a "surrogate" series is generated. Because of the 



Wiener-Kinchine theorem 



401 ]. the surrogate record has the same linear correlations as the 



original, as well as the distribution. The only difference is that the original record has 
the non-linear correlations (if they exist) but the surrogate does not have any non-linear 
features. 

In Fig. [1] we also plot the cumulative distribution for the surrogate with the same three 
thresholds as the original. Since the surrogate records lost the non-linear correlations, they 
are similar to each other, we only show results for GE's surrogate. It is seen that the collapse 
of the surrogate for different q values is significantly better than that of the original and the 
deviation tendency with the threshold in the original records disappears. This indicates that 
the scaling deviations in the original are due to the non-linear correlations in the volatility. 
To further test this hypothesis, we analyze the moments [L m = ((t ■/ '{V)) m ) 1//m in Sec. (IV) 
and show similar and consistent deviations from scaling. We also compare our results to the 
stretched exponential distribution (dashed lines). This curve is very close to the empirical 
results, in particular for the surrogate records which contain only the linear correlations. 



5 



This suggests that PDF of return intervals is well approximated by a stretched exponential. 



IV. THE MOMENTS OF SCALED INTERVALS 



The cumulative distribution shows clear systematic trend with q, which is difficult to see 



from the PDF directly 



22 



23 



24. 



25 



26. 



271 ] . To further analyze the systematic tendency in 



the distribution, we calculate the moments fi m averaged over a stock dataset as a function 
of (r), where a mean interval (r) corresponds to a threshold q and therefore characterizes a 
return interval series. We study moments for a wide range of (r), from 3 minutes (to avoid 
the artificial effects due to discreteness close to r = 1) to thousands minutes (few trading 
days or even a week). Assuming a single scaling function for the PDF P g (r), Eq. (JTJ, it 
follows 

/*oo 1 /*oo 

A*m = ((r/(r)D 1/m = { / (r/(r}) m - T -f(r/(r})dr} 1 / m = { x m f(x)dx}^ m , (7) 

Jo v) Jo 

which only depends on m and on the form of the scaling function f(x) but independent of 
(r). Thus, if \i m depends on (r), it suggests deviation from the assumption of scaling. 



A. Moments vs. Mean Interval (r) 

First we examine the relation between the moments fi m and the mean interval (r). Fig. [2] 
shows four representative moments m = 0.25, 0.5, 2 and 4 for stock C, GE, KO and XOM. 
Ignoring small fluctuations, which is usually due to limited size data, all moments \x m for the 
original records deviate significantly from a horizontal line, which is expected for a perfect 
scaling of the PDF. They depend on (r) and show some systematic tendency. For m > 1, 
moments have similar convex structure, first \x m increases with (r) and then decreases, where 
the crossover starts earlier for larger m. For m < 1, moments also show similar tendency 
but in the opposite direction compared to m > 1. These deviations from scaling in fi m 
are consistent with the deviations seen in the cumulative distributions shown in Fig. [IJ 
Moments of large m (m > 1) represent large r/(r) in the PDF and they initially (for 
(t) < 100) increase with (r), while moments of small m {m < 1) represent small rj (r) and 
they initially (for (r) < 100) decrease with (r). 

To further test if the systematic deviations are not due to finite size effects and dis- 
creteness, we also examine moments for the surrogate records which are more flat for most 
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range, as shown in Fig. [21 For the same order m, the moment of the surrogate obviously 
differs from that of the original, especially in the medium range of (r) (10 < (r) < 100). 
This discrepancy suggests that the non-linear correlations exist in the original volatility 
and accounts for the scaling deviations. Nevertheless, all moments of surrogate show small 
curvature from a perfect straight line at both short and long (r), which are much weaker 
compared to the original records. The weak curvature suggests that some additional effects, 
not related to the non-linear correlations, affect the moments. For small (r), the resolution 
discrete limit seems to have some influence on the moments. We will discuss this effect in 
section C. For large (r), the moments are gradually approaching the horizontal line and are 
more fluctuating, the effect seems to be related to limited size of the record. This effect will 
be discussed in section B. 

B. Multiscaling 

For the original volatility records, the systematic tendency in the distribution of r and the 
moments implies that the return intervals may have multiscaling features. To avoid as much 
as possible the effect of discreteness and finite size, we calculate the moments only for some 
medium range of (r) where the effects are small. Since there is no non-linear correlations in 
the surrogate records, the curvature in their moments is only due to the additional effects, 
we use the surrogate curve as our reference. For small (r), the increasing (decreasing) range 
for m > 1 (m < 1) almost ends at (r) = 10 minutes. For large (r), the curves start to 
decrease (increase) from different positions, but at (r) = 100, all curves do not or just 
start to decrease (increase). Thus we choose to study \i m in the region, 10 < (r) < 100, 
represented by the shadow areas in Fig. [2j In this range, we find a clear trend for the 
original records while the surrogate is almost horizontal. To quantify the tendency, we fit 
the moments with a power-law, 

~ (r) a - (8) 

If the distribution of rj (r) follows a scaling law, the exponent a should be some value very 
close to 0. If a is significantly different from 0, it suggests multiscaling. 

To examine the multiscaling behavior for the whole market, we calculate a for all 500 
stocks of S&P 500 constituents and plot the histogram for m = 0.25 to 2. Fig. [3] shows 
that each histogram has a narrow distribution, which suggests that a are similar for the 500 
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stocks. For the original records, almost all a significantly differ from 0, thus the moments 
clearly depend on the mean interval. Moreover, the mean value of a shifts with order m 
from (a) ~ —0.2 for m = 0.25 to (a) ~ 0.1 for m = 2 which means the dependence varies 
with the order m. This behavior suggests multiscaling in the return intervals distribution. 
Indeed, histograms for the surrogate records are more centered around values close to a = 0. 
The uniscaling behavior for the surrogate suggests that the non-linear correlations in the 
volatility are responsible for the multiscaling behavior in the original. 

To remove fluctuations and show the tendency clearer, we plot the dependence of (a) 
on m, where (a) is the average a over all 500 stocks. In Fig. H] we show this relation for 
a wide range of m, 0.1 < m < 10, and the plot shows two different behaviors. For small 
m (roughly m < 2), (a) for the original records clearly deviates from and demonstrates 
the multiscaling behavior, while (a) for the surrogate is closer to 0. For large m (m > 2), 
the two curves have similar decreasing trend. Since large r dominates high order moments, 
this similarity may be due to finite size effects. To test the finite size effects we simulate 
surrogate return intervals by assuming a stretched exponential distribution i.i.d. process 
with 3 sizes (number of all r in the series), 2 x 10 6 , 2 x 10 5 (the size of the empirical dataset) 
and 2 x 10 4 . Without loss of generality, we choose 7 = 0.3, which is the correlation exponent 
for GE of q = 2. To be consistent with the 500 stocks, we perform 500 realizations and 
plot their average exponent (a). As shown in the inset of Fig. HI all (a) show a similar 
decreasing trend as that of the empirical curves. However, it is seen that the trend starts 
earlier for smaller size, and thus the size limit has a strong influence on high order moments. 
Fig. H] also shows the error bars for the two records, which is the standard deviation of 500 
a values. Note that the error bars for the volatility records do not overlap those of their 
corresponding surrogate, indicating the significance of our results. 

C. Discreteness Effect 

For small (r) ((r) < 10), the behavior of function of (r) was attributed to the 

discreteness. Here we examine this effect. Due to the limits in recording, we can not have 
a continuous but discrete record. In our study the volatility is recorded in 1 minute. The 
relative errors in moments will be considerable large for small (r) close to 1 minute. By 
starting from (r) = 3, we only partially avoid the discreteness in the moments. To test the 
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discreteness effects, since we can not increase the resolution, we reduce it and compare the 
moments fi m with 3 resolutions, 1 minute, 5 minutes and 10 minutes. Fig. 0(a) shows this 
comparison for GE for m = 0.5 and m = 2. The three resolutions have a similar trend, 
showing that the curves become flatter for the higher resolution. For other stocks, we find 
similar behavior. This systematic tendency suggests that the recording limit (1 minute) 
strongly affects the moments at short (r). To reduce it, we should raise the recording 
precision or study the moments of larger (r). 

To further test this result, we simulate artificial return intervals with an i.i.d. process 
from stretched exponential distribution with 7 = 0.3, same as the empirical 7 for GE of 
q = 2. We examine moments of m = 0.5 and m = 2 with the same three resolutions (1, 5 
and 10 time units) as in the empirical test done above. The simulated size is 200 thousands 
points for each trial and we use the average over 100 trials for each resolution. Fig. [5(b) 
shows curves similar to that of empirical (Fig. [5(a)). For the higher resolution, the curve 
is closer to the horizontal line and finally may reach the line when we raise the resolution 
high enough. To show this, we also simulate continuous return intervals and find constant 
moments, as expected (Fig. [5(b)). Therefore, the discreteness effects can be overcome if the 
resolution is improved enough. Note that for the empirical data, we expect the moment 
not to be constant for small (r) even if we have a much better resolution, since the return 
intervals has the multiscaling behavior, as shown for larger (r) in the range 10 < (r) < 100 
which is not affected by discreteness. 

D. Moments vs. Order m 

The moments \i m have systematic dependence on m, as seen in Figs. [2] and [3] where the 
moments are plotted as the function of (r). It is of interest to explore the relation between 
the moments and m directly. For a fixed (r), representing a given threshold q one can 
study, the return intervals and their moments of various orders which exhibit information 
on different scales of r. Moments of large m represent large r and vice versa. If rj (r) follows 
a single distribution without corrections due to effects such as discreteness and finite size, 
curves of \x m vs. m for different (r) should collapse to a single one, which only depends on 
the scaling function f(x) from Eq. ([7j). In Fig. [6] we plot /i m vs. m for both the original and 
surrogate records. We plot \x m for m between 0.1 and 10 for three (r) values: 10, 80 and 
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400 minutes. For the original (Fig. E^a)), there is substantial deviations from a single curve. 
This supports our suggestion that the return intervals has multiscaling behavior. Moments 
for the surrogate (Fig. E^b)) converge to a single curve for m < 2 but become diverse for 
high orders, which agrees with the strong influence of the finite size effects. As a reference, 
we also plot the analytical moments (Fig. M^c)) from the stretched exponential distribution. 
Substituting Eq. ([2]) into Eq. ([7]), we obtain 

^ m ~a { r(i/ T ) 1 ■ [J) 

Fig. E](c) shows analytical curves for various correlation exponent 7. 
V. DISCUSSION 

We study the scaling properties of the distribution of the volatility return intervals for 
all S&P 500 constituents. We find small but systematic deviations from scaling assumption 
with the threshold q in the cumulative distribution. Compared to the good collapse for 
the surrogate records where non-linearities are removed, this suggests that the origin of 
this trend is due to non-linear correlations in the original volatility. Moreover, we find 
similar systematic deviations for the moments fi m , which are also attributed to the non- 
linear correlations in the volatility. We distinguish these deviations from the deviations due 
to the discreteness for small (r) and finite size effect for large (r). Further, we explore the 
dependence of the moment /i m on its order m. When compare to surrogate records and 
to analytical curves, the results support the multiscaling hypothesis of the return intervals. 
Thus, the scaling assumption in the return interval distributions although it is a good 
approximation can not be exact. Also, the stretched exponential form of the scaling function 
can only be an approximation. Recently Eisler et al. 32|, |33| exhibits that the distribution 
of intertrade times has similar multiscaling behavior and the market activity depends on 
the company capitalization. It would be interesting to connect the intertrade times with the 
return intervals and test size dependence in the return intervals. 
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x/<x> 

FIG. 1: (Color online) Cumulative distribution of the scaled intervals t/(t) for stock C, GE, KO, 
XOM and GE's surrogate. Symbols are for three thresholds q = 2 (circles), 4 (squares) and 6 
(triangles) respectively. As two examples, we fit the cumulative distribution of stock C (q = 2) and 
the surrogate (q = 2) to a stretched exponential distribution (Eq. fl2J)) with exponent 7 = 0.25 and 
7 = 0.50 correspondingly. Except that for stock C, all symbols and curves are vertically shifted 
for better visibility. For the surrogate records, symbols collapse almost perfectly to one curve. 
However, all original data exhibit similar systematic deviations from the scaling. This suggests the 
non- linear correlations in the volatility series affects the scaling in its return intervals. 
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FIG. 2: (Color online) Moment \x m for the scaled intervals of stock C, GE, KO and XOM. We 
show results from the original volatility series (filled symbols) and their surrogate (empty symbols). 
For each case, four moments, m = 0.25 (circles), 0.5 (squares), 2 (diamonds) and 4 (triangles) 
are demonstrated. Symbols for the original are clearly away from the horizontal line and their 
deviations are much larger than that for surrogate, therefore the non-linear correlations in the 
original are related to those deviations. To avoid effects from the resolution and size limit, we 
choose the shadow area, 10 < (r) < 100, to study the multiscaling behavior. 
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FIG. 3: (Color online) Distribution of multiscaling exponent a for S&P 500 constituents. The 
exponent a is obtained from the power-law fitting for moments in the medium range 10 < (r) < 100. 
(a) Histogram of a for the original volatility and (b) for surrogate. The distributions have a 
systematic shift with m in (a) while all of them almost collapse in (b). This suggests the multiscaling 
behavior in the return intervals in the original records. The significant discrepancy between the 
original and the surrogate records manifests that non-linear correlations form the original volatility 
accounts for the multiscaling behavior. 
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FIG. 4: (Color online) Dependence of average multiscaling exponent (a) on order m. The average 
(a) was taken over the 500 a of the S&P 500 constituents. Results for the original (circles) and 
surrogate records (squares) are displayed. For large m, the two curves have the similar tendency 
which is attributed to the finite size effects. For small m, the two curves are significantly different 
which supports the multiscaling in the return intervals, due to the non-linear correlations in the 
original volatility. The inset demonstrates {a) averaged over 500 stretched exponential distributed 
i.i.d. simulations with 7 = 0.3. Three sizes, 2 x 10 6 , 2 x 10 5 and 2 x 10 4 are displayed, which 
clearly shows the finite size effect. 
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FIG. 5: (Color online) Discreteness effect in the moments, (a) Moments for stock GE with three 
resolutions, 1 (circles), 5 (squares) and 10 minutes (triangles). Filled symbols are for m = 0.5 while 
empty symbols are for m = 2. (b) Moments of artificial records averaged over 100 simulations. 
For each trial, we simulate the return intervals which follows a stretched exponential distribution 
with 7 = 0.3 and the length of 200 thousands points. Symbols are similar to that in (a). Three 
resolutions, 1, 5 and 10 time units are displayed. To show the disappearing of the discreteness, we 
also plot simulation results for continuous return intervals (triangles), which exhibits a constant 
moment, independent of (r). 
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FIG. 6: (Color online) Dependence of the moment \i m on the order m. (a) for the return intervals 
of the original volatility for stock GE. (b) for its surrogate. Three mean interval (r) = 10 (circles), 
80 (squares) and 400 minutes (triangles) are demonstrated in (a) and (b). (c) Analytical moments 
from stretched exponential distributions, 7 = 0.25, 0.50, 0.75 and 1 (Poisson distribution) taken 
from Eq. Q. For large m, both (a) and (b) show discrepancies which is related to the finite size 
effect. For small m, the difference between (a) and (b) is due to the non-linear correlations in the 
original volatility. Compared to the analytical curves in (c), the return intervals for the original 
records shows multiscaling behavior. 
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